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We report on the extraction of the target mass contributions to the unpolarized proton structure 
functions by applying an unfolding procedure to the available world data from charged lepton scat- 
tering. The method employed is complementary to recent and future parton distribution function 
fits including target mass contributions and the results obtained can be utilized to further study per- 
turbative QCD at large Bjorken x and small Q 2 . Global fits are performed to both the available F% 
and the separated F£ (Ff) data which yield excellent descriptions and provides parameterizations 
of the extracted target mass contributions. 

PACS numbers: 



I. INTRODUCTION 

Since the inception of QCD as the theory of the strong interaction, perturbative QCD (pQDC) has been spectacularly 
successful in explaining a host of scattering and annihilation processes involving hadrons. Among these is the Q 2 
dependence of the proton F 2 structure function measured in charged lepton deeply inelastic scattering (DIS). However, 
the separation of the pure leading twist pQCD Q 2 dependence from possible higher twist (HT) and so called target 
mass contributions (TMCs) is required for any quantitative study of pQCD at moderate to large values of Bjorken 
x and low Q 2 . The large body of existing data on F% covers a significant range in Q 2 and x, and is a critical input 
for pQCD based fits which attempt to extract the parton distribution functions (PDFs). While the contributions 
from TMCs are now being included in some PDF fits, cf. e.g. Historically the data which were expected to 

have significant TM or HT contributions were excluded in standard analyses from the fits via kinematic cuts on the 
invariant hadron mass, W. Typically, this has removed nearly all the precision data from SLAC at moderate to large 
x from the fits, significantly reducing the Q 2 lever arm on the remaining data, and resulting in the rather large error 
bands quoted from these fits for x > 0.5. 

In this paper we present a fit to the world charged lepton scattering data on the proton F% and Fl (F±) structure 
functions with the aim of separating the TMCs from the massless limit structure functions via an unfolding procedure. 
The method employed is complementary to PDF fits including the twist-2 target mass contributions but has the 
advantage that it does not rely on a finite expansion in the strong coupling constant a s . This avoids the issue of 
disentangling higher order QCD corrections from the TMCs. The extracted distributions still contain the higher twist 
terms. 

II. DATA SETS 



The world charged lepton scattering data sets utilized for the F 2 fit are listed in Table U and include data from 
the BCDMS [7j, and NMC [8[ collaborations, the reanalyzed SLAC data from Whitlow et al. [9J, and the e + and e~ 
data from the HI [l(| [H| and ZEUS [H, [l3[ collaborations; the latter two help to constrain the fit at both low x and 
large Q 2 . To help constrain the low Q 2 region, recent data from Jefferson Lab Hall C experiments E99-118 [IH and 
E94-110 [Hi have been utilized. 

These experiments have additionally performed limited longitudinal and transverse (L/T) separations to extract 
the information on the transverse structure function, F\, the longitudinal structure function, Fl, and the ratio 
R = Fl/2xFx. For the Fl fit we include only the data extracted via L/T separations in order to minimize the 
correlation between the F 2 and Fl data sets. For the JLab experiments, only the L/T separated values for F 2 were 
utilized. It is unfortunate that more experiments have not made readily available the complete set of separated 
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structure functions extracted from L/T separations. Instead, what can generally be found is a table of the extracted 
R = (Jl/^t = Ft/2xFi values, but without the consistently extracted values for F 2 (or Fx). In these cases the 'data' 
on Ft have been constructed utilizing the available data on R and a parameterization of the F 2 structure function. 
The longitudinal structure function was then constructed from 

r 2 F 2 

with r 2 = = 1 + For this purpose the phenomenological ALLM fit [l6[ was utilized, which reproduces the 

world F 2 data to better than 3% and we include this as an additional uncertainty when propagating the uncertainties 
on Fl (Fi), which are ultimately dominated by the uncertainties on the R data. The data sets utilized for the Ft, fit 
are listed in Table HU 

Before discussing the procedure, we feel it is worth commenting further on the need for uncorrelated data sets for 
Ft and F 2 . This is also important for extractions of PDFs utilizing pQCD fits, since Ft provides direct information on 
the gluon content that is only present in F 2 through the shallow logarithmic scaling violations. In principle, the most 
information on the structure can be obtained from the full data set on the reduced cross section, which is proportional 
to r 2 F 2 — (1 — e)Ft(x, Q 2 ), with e a kinematic factor describing the photon polarization. Most of the DIS data sets 
extract F 2 utilizing all the cross section measurements, after first determining R by a linear fit in e to the reduced 
cross section measurements at different e values at fixed x and Q 2 (the Rosenbluth separation technique [l7j). 

III. PROCEDURE 

We follow the formalism for including TMCs to the massless limit structure function developed by Georgi and 
Politzer |6| within the operator product expansion (OPE). This approach was extended by Kretzer and Reno [18| to 
include the full set of electro- weak structure functions using modern formalism. For a recent review of TMCs see [19j ■ 

At fixed Q 2 the F 2 structure function, including TM effects, can be written 



f 

Here, F 2 °\x, Q 2 ) is the massless limit structure function which can be calculated in pQCD utilizing the PDFs 
determined at some scale, £ = 2x/(l + r) is the Nachtmann scaling variable, and both x' and x" are dummy 
integration variables. The corresponding TMC equation for Fl is 

fjw) = ^V 2) + ^4 /V^£^ + g£< /V /V ^ Q2) - (3) 



x' 2 Q i J ( J T , X " 2 



e Q 

The transverse structure function including the TM contributions, can then be constructed from 



and in the massless limit (r — > 0) from 



rpTM jpTM 

2xF TM = ^2 H_ (4) 



2xF 1 (0) =F 2 (0) -Ff ) . (5) 



We note that these equations are valid to all orders in pQCD at twist 2. The target mass corrections for for higher 
twist operators in the unpolarized case are yet unknown. Moreover, dynamical higher twist contributions, starting 
with twist 4, cannot be described with one scaling variable only, neither do their potential target mass corrections, 
because the main variable x is supplemented by further other invariants Xi. Their number grows with the level of 
twist included. This implies the necessity to study different multi-parton correlations, which cannot be determined 
by inclusive measurements as that of the deep-inelastic structure functions only. Because of this the present analysis 
remains phenomenological applying the twist 2 target mass corrections to the structure functions. [30j | 

While it is F 2 M (x, Q 2 ) that is measured in scattering experiments, it is the Q 2 evolution of F^ that is described 
by pQCD. To extract F^ from data requires the inversion of Eq. ([2]). The procedure for extracting the TM corrected 
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F 2 structure function is to parameterize F$ \x 7 Q 2 ) and then to minimize the x 2 difference between the left-hand 
side of Eq © and the F 2 data: 

/ J?TM jpi\ 2 

where F 2 is the value of the i th data point and SF 2 it's uncorrelated uncertainty. 



TABLE I: Data sets used in F 2 fit. 



Data Set 


Qtlin 




QlAax 


%max 


# Data Points 




(GeV 2 ) 




(GeV 2 ) 






BCDMS [2 


7.5 


0.070 


230 


0.75 


178 


NMC [8] 


0.75 


0.0045 


65 


0.50 


158 


SLAC (Whitlow [9]) 


0.58 


0.063 


30 


0.90 


661 


HI [10,^1] 


1.5 


3.0xl0 -5 


5000 


0.32 


286 


Zeus [12, ,13] 


3.5 


6.3 x 10" 5 


5000 


0.20 


228 


E99-118 [14] 


0.273 


0.077 


1.67 


0.320 


9 


E94-110 [15] 


0.7 


0.19 


3.5 


0.57 


12 



TABLE II: Data sets used in F L (Fi) fit. 



Data Set 


(GeV 2 ) 




Qlvlax 

(GeV 2 ) 


%max 


# Data Points 


BCDMS [7] 


15 


0.07 


50 


0.65 


10 


EMC [21] 


15 


0.041 


90 


0.369 


28 


NMC [8] 


1.31 


0.0045 


20.6 


0.11 


10 


SLAC (Whitlow [22]) 


0.63 


0.1 


20 


0.86 


90 


SLAC (E140x [23]) 


0.5 


0.1 


3.6 


0.50 


4 


HI [24] 


1.5 


3.0 x 10" 5 


45 


0.0015 


13 


ZEUS [25] 


24 


6.7 x 10" 4 


110 


0.0049 


18 


E99-118 [26] 


0.273 


0.077 


1.67 


0.320 


8 


E94-110 [15] 


0.7 


0.19 


3.5 


0.57 


12 



The fit form utilized for the x dependence of both F 2 and F£ was 

F^l{x) = Ax B (l-xf(l + DVx + Ex), (7) 
with the Q 2 dependence of A, B , C, D, and E parameterized as 

A(Q 2 ) = A 1 + A 2 e-Q 2 / As + Ai log(0.3 2 + Q 2 ). (8) 

These forms are empirical, but were found to give a good description of the data. The fitting was performed utilizing 
MINUIT [27J and utilized an optimized FORTRAN routine DAIND 29] to perform the single and double integrations in 
Eqs. © and ©. 



IV. RESULTS AND CONCLUSIONS 



The fit results for both F 2 and Fl are quite good, with x 2 P er degree of freedom (x 2 /DoF) of 1.01 and 0.93 for F 2 
and Fl, respectively. The resulting fit parameters are listed in Table ILTTl for F 2 and Table ITVl for F^. Comparisons 
of the fit results to the F 2 data are shown in Figure [T] for the 6 different Q 2 values listed on the panels from 1 to 
100 GeV 2 . Data within a range of Q 2 about the central value have been centered to the common Q 2 value using 
the fit. A comparison over the full x range fit is shown (Left), as well as for large x (Right). Similar comparisons 
for Fl are shown in Figure [5] Of particular note is the lack of quality data for Fl in the range 0.005 < x < 0.3 for 
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TABLE III: F 2 best fit parameters. 



n — 


1 


2 


3 


4 


A' 2 


0.2145 











Bl 


-2.1540 


.0808 


5.6730 


-0.0319 


c„ 


2.7830 


-0.5013 


7.104 


0.1857 


D 2 


-0.8997 


-3.1092 


0.3820 


-0.1275 


E n 


9.0441 












TABLE IV: Fl best fit parameters. Parameter (indicated by *) has been fixed, as the goodness of the fit was found to be 
only minimally sensitive to it's value. 



n = 


1 


2 


3 


4 




0.0078 


0.0 


0.0 


0.0 


B n 


-2.0650 


0.0 


0.0 


-4.8846 


r L 


1.0* 


0.0 


0.0 


0.0 


D L n 


0.0 


-0.0030 0.0 


0.6040 




0.0 


0.0 


0.0 


0.0 



Q 2 > 3 GeV 2 . Such data is critical for constraining the fit at low x and would be of particular value in constraining 
the gluon distribution which is a dominant contribution to F^. 

We expect that the fit will be of interest to the broader community in providing a good description of the data on 
the full set of unpolarized structure functions in charged lepton scattering from the proton, F2, Fl, F\, and 7?, over 
the kinematic range of x > 3 x 10~ 5 and 0.1 < Q 2 GeV 2 , in addition to providing a separation of the TMCs from the 
massless limit structure functions. 

The separated structure functions can be used for further pQCD studies, such as HT analyses of both the structure 
functions and moments. The current analysis procedure is complementary to PDF fits which include TMCs in that 
the unfolding of the TMCs does not depend on a finite order expansion in as. Comparison of the results from the 
two procedures could provide an additional handle on isolating HT contributions from higher-order QCD corrections 
and making for a 'cleaner' evaluation of the former. 

Additionally, the fit will be of use in future studies of quark-hadron duality. In particular, determining how well 
the averaged resonance region structure functions at fixed x obey the Q 2 dependence dictated from DIS. 

FORTRAN computer code embodying the fit described in this article is available by email request from the authors, 
including the full covariance matrices. 
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FIG. 1: (Color online) Comparison of the fit results (dashed curve) to the proton data for 6 different Q 2 values (left 

panel). The low x data from HI and ZEUS are beyond the vertical scale. Also shown is the massless limit structure function 
parameterization from the fit (solid curve). A zoom view for the region x > 0.5 is also shown for 5 different Q 2 values (right 
panel) . 
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FIG. 2: (Color online) Comparison of the F£ fit results (solid curve) to the proton data versus x for 5 different Q 2 values (left 
panel). The low x data from HI and ZEUS are beyond the vertical scale. Also shown is the massless limit structure function 
parameterization from the fit (dot-dashed curve). A zoom view for the region x > 0.5 is also shown (right panel). All data 
within a range of ±40% (Q 2 = 1), ±33% (Q 2 = 3), ±15% {Q 2 = 5, 7, 10, and 12), and ±20% (Q 2 = 20 and 40) have been 
bin-centered to the central Q 2 using the Q 2 dependence from the fit. The uncertainty band stemming from the uncorrelated 
uncertainties and calculated from the fit covariance matrix is indicated by the dashed curve. 
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